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Tight-binding molecular dynamics (TBMD) is used to 
study the structural and electronic properties of silver clus- 
ters. The ground state structures of Ag clusters up to 21 
atoms are optimized via TBMD combined with genetic algo- 
rithm (GA). The detailed comparison with ab initio results 
on small Ag„ clusters (n = 3 — 9) proves the validity of the 
tight-bind model. The clusters are found to undergo a tran- 
sition from "electronic order" to "atomic order" at n = 10. 
This is due to s-d mixing at such size. The size dependence 
of electronic properties such as density of states (DOS), s- 
d band separation, HOMO-LUMO gap, and ionization po- 
tentials are discussed. Magic number behavior at Ag2, Agg, 
Agi4, Agis, Ag2o is obtained, in agreement with the predic- 
tion of electronic ellipsoid shell model. It is suggested that 
both the electronic and geometrical shell exist in the coinage 
metal clusters and they play a significant role in determining 
cluster properties. 
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I. INTRODUCTION 

The structural and electronic properties of metal clus- 
ters are currently a field of intensive research both the- 
oretically and experimentallyErEl. The basic theoretical 
concept in the electronic structure of alkali metal clusters 
is the shell mfMijl based on jellium sphere (or ellipsoid) 
approximatiorua, which has successfully interpreted the 
magic number effect in Na„ and K„ clusters {n =2, 8, 20, 
40, • • ■). As compared to alkali-metal clusters, the appli- 
cation of electronic shell model to coinage-metal clusters 
(Cu„, Ag„, Au„) is somewhat questionable due to the 
existence of inner d electrons. Among the noble metal 
clusters, Ag„ is expected to exhibit the largest similarity 
to the alkali metal clusters as the 4c? orbitals in Ag atom 
are low- lying and act almost like innershell core orbitals. 

Experimental studies on silver clusters include mass- 
spectra i ,|-| ionization potentials (IPs)bIj, rfihotoelectron 
spectraOia, electrcp_spin resonance (ESR)ll3, optical res- 
onance absorptioniJ, etc. In general, most of the cluster 
properties resemble the predictions of shell model within 
one s electron picture. But there are still some experi- 
mental evidences for Ag„ that are different from alkali- 
metal clusters and cannot be understood via the shell 
model of s electrons. For instance, the Mie resonance 
peak of silver clusters exhibit blue shift with decreas- 
ing cluster radiusO, while red shift of Mie resonance fre- 



quency is found for alkali-metal clusters. A direct com- 
parison of photoelectron spectra between Ag„ and Na„ 
display signifijCant difference that might be attributed to 
the d orbitalsB. Therefore, it is important to clarify the 
contribution of id electrons and s — d interaction in the 
silver clusters and the geometrical effect on the electronic 
properties of the clusters. 

Besides shell model, the metal clusters have been inves- 
tigated by accurate quantuni, chemical approaches based 
on molecular orbital theoryu. However, such ab initio 
calculations on coinage metal clustecs |a£e quite time con- 
suming and limited in small sizeE3iij. Among those 
works, the most detailed and comprehensive study of 
small neutral silver clusters (Ag2 to Agg) has been per- 
formed via configuration interaction (CI) noathod with 
relativistic effective core potential (RECP)l3. On the 
other hand, the electronic structures of several larger 
silver clusters has beep-, approximately described by a 
modified Hiickel modelES. However, all these studies are 
carried out for limited number of structural candidates 
with symmetry constrain. An unbiased minimization for 
the cluster ground state structure incorporated with elec- 
tronic structure calculations would be much more infor- 
mative for understanding the interplay between geomet- 
rical and electronic structure and testing the validity of 
electronic shell model. 

Up to now, the most reliable and accurate procedure 
in dynamic searching of cluster equilibrium geometries is 
provided by Cas-Parrinello (CP) method with simulated 
annealing (SA)tZI. But such kind of ab initio simulation 
is also limited in small size (about n < 10) in a truly 
global optimization because of the rapidly increase in 
computational expense with cluster size. Among coinage 
metal clusters, the CP method has been emplojied to 
study small Cu„ clusters {n = 2 — 10) recentlyll3, but 
the corresponding investigation on silver and gold clus- 
ters is not available so far. In recent years, tight-binding 
molecular dynamics has been developed as an alternative 
to CP method in atomistic simulation for larger scaled 
systemst3. As compared to ab initio methods, the pa- 
rameterized tight-binding (TB) Hamiltonian reduces the 
computational cost dramatically. It is more reliable than 
classical simulation based on empirical potential since the 
atomic motion is directly coupled with electronic struc- 
ture calculation at each steps. For transition- metal clus- 
ters, M.Menon and co-workers have proposed a minimal 
parameter tight-binding_.acheme and used it to study 
nickel and iron clusteraHj^El. In this work, we shall in- 
troduce a similar TB model for silver. By using the 
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TBMD as aJotal minimization approach, genetic algo- 
rithm (GA)E3~E3 is employed to search the global min- 
imal structure of Ag„ clusters up to 21 atoms. The 
size dependence of relative stabilities, density of states 
(DOS), HOMO-LUMO gaps, and ionization potentials 
of the clusters are calculated and compared with avail- 
able experimental results. The magic effect of electronic 
shell and the interplay between geometrical and elec- 
tronic structure in silver clusters are discussed. 



II. THEORETICAL METHOD 

In the minimal pararrLcter tight-binding scheme pro- 
posed by M.Menon et altB, the total binding energy Eb 
of transition-metal atoms is written as a sum: 



Eb — Eel + Er 



(1) 



Ef,i is the electronic band structure energy defined as the 
sum of one-electron energies for the occupied states 



occ 

Eel = ^ fife 
fc 



(2) 



where energy eigenvalues Cfc can be obtained by solving 
orthogonal 9n x 9n TB Hamiltonian including Ad, 5s and 
5p electrons. The repulsive energy Erep is described by 
a pair potential function x(r) of exponential form: 



-4:Oi.{r — dQ) 



(3) 



j>i 



where is the separation between atom i and i. dn = 



2.89A is the bond length for the fee bulk silveio, a is 
taken to be one-half of l/dg according to Ref.20. 

In order to reproduce the cohesive energies of small 
clusters through bulk TB hopping parameters, it is still 
necessary to introduce a bond-counting term Ebond- 



Ebond = ~N[ainb/N) + b] 



(4) 



Here the number of bonds Ub are evaluated by summing 
over all bonds according to cutoff distance Rc 



Ub 



[exp(- 



Rr 



(5) 



Vxx',id) = Vxx'Mi^) 



r+2 



(6) 



The parameter t — for s-s, s-p, p-p interactions, r = 
3/2 for s-d and p-d interaction, r = 3 for d-d interaction. 

In present, we have adopted the Slater-Koster hopping 
integrals Vxx'fj.ido) and the on-site orbital energy from 
the values fitted to first perinciple APW band structure 
calculation of bulk silverEZl. Furthermore, to transfer the 
on-site orbital energies levels from bulk calculation to 
atomic limit, a constant energy shift Ae=-15.88 eV is 
applied on the on-site energies from Ref.27. Such shift 
in on-site levels make the theoretical ionization potential 
of Ag„ clusters quantitatively comparable to experimen- 
tal values. The repulsive potential parameter xo is fitted 
for experimental bond length 2.48 A of silver dimerEa. 
The bond-counting terms a, b are chosen to reproduce 
the Jib^initio binding energy for small clusters Ag3 , Ag4 , 
Ag5E3a'E£l. All the parameters used in our calculation 
are listed in Table I. These empirical parameters can de- 
scribe both bulk phase and dimer of silver with an accept- 
able accuracy. The cohesive energy 2.75 eV, equilibrium 
interatomic distance 2.88 A of fee solid silver from TB 
model are close to the corresponding experimental value 
2.95 eV and 2.89 A respectivelyEa. The vibrational fre- 
quency and binding energy calculated for silver dimer at 
equilibrium distance is 125.5 cm~^ and 1.25 eV, in rea- 
sonable agreement with experimental data of 192.4 cm~^ 
and 1.66 cVEI. 



Table I. 
work. 



Parameters of TB model for silver used in this 



-6.452 eV 



-0.447 eV 



-14.213 eV 



-14.247 eV 



VjjptT VppTT VsdcT 



-0.895 eV 1.331 eV 2.14317 eV 0.088 eV -0.423 eV 

Vpdcr VpdTT Vdda Vddn VddS 

-0.531 eV 0.207 eV -0.429 eV 0.239 eV -0.046 eV 



do a xo o. b Rc A 

2.89 A 0.692A"1 0.58 eV -0.16 eV 0.59 eV 3.5 A 0.1 A 



It should be noted that only the first two terms Eei and 
Erep in Eq.(l) contribute to the interatomic forces in 
TBMD simulation, while the Ebond term is added after 
the relaxation has been achieved. However, for metal 
clusters, this correction term is significant|iB distinguish- 
ing various isomers at a given cluster sizeE2l. 

The 9?! X 9?! TB Hamiltonian matrix is constructed 
with Slater-Koster scheme, while the distance scaling of 
hoppingjintegrals Vxx'iJ, is taken as the Harrison's univer- 
sal formEZi: 
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The determination of lowest energy structuc^pOf clus- 
ters is performed by a genetic algorithm (GA)c3i£3. The 
essential idea of this novel global optimization strategy 
is to mimic the Darwinian biological evolution process 
in which only the fittest candidate can survive. Some 
pioneering works have demonstrated the impressive ef- 
ficiency of GA in searching the global minima of clus- 
ters as compared to standard simulated annealing. At 
beginning, we generate a number Np of initial configura- 
tions by random {Np = 4 — 16, depending upon cluster 
size). Any two candidates in this population can be cho- 
sen as parents to generate a child cluster through mat- 
ing process. In the end of mating procedure, mutation 
operation is allowed to apply on the geometrical struc- 
ture of child cluster with 30% possibility. The child clus- 
ter from each generation is relaxed by TBMD quenching 
of 300-500 MD steps with At 5/s. Then the locally 
minimized child is selected to replace its parent in the 
population if it has different geometry but lower binding 
energy. Typically, 100—150 GA iterations is sufficient to 
ensure a truly global search up to n = 21. The detailed 
descripticpn|arbout the practical operation of GA is given 
elsewhereE3e3. 



III. STRUCTURES AND STABILITIES OF 
SLIVER CLUSTERS 

A. Structures of small Agn with n < 9 

By using the combined GA-TBMD strategy, we have 
explored the global minimal structures of Ag„ clusters up 
to 21 atoms. These ground state structures are shown 
in Fig.l (4 < n < 9) and Fig.2 (10 < n < 21). In 
Table II, the structural parameters (bond lengths and 
bond angles), binding energies and ionization potentials 
of the small Ag„ clusters (n = 3 — 9) in ground state and 
some metastable structures arCpi^Dnupared with accurate 
quantum chemistry calculationaiJilj. The lowest energy 
structures found for most of the clusters coincide well 
with mteysive CI calculations in Ref.l4 and the other 
workal3£3^E£l. The calculated cluster properties of ground 
state and metastable isomers agree with ah initio results 
quantitatively. As shown in Table I, the TB bond lengths 
are typically within 0.05—0.15 A according to the ah 
initio values. The average deviation of binding energy per 
atom and iopi^ation potentials from this work to ah initio 
calculationsliJ is about 0.13 eV and 0.30 eV respectively. 




FIG. 1. Lowest-energy structures for Ag„ (n = 4 — 9) clus- 
ters. 

For silver trimer, isosceles triangle is about 0.05 eV 
lower in energy than equilateral triangle and 0.84 ey. 
lower than the linear isomer. ESR experiments on AgsEI 
supports the isosceles triangle structure with C2v symme- 
try. In the case of Ag4, planar rhombus is lower in energy 
than a relaxed tetrahedron by Ai?=0.31 eV although 
the tetrahedron has higher coordination and has been 
predicted as most st|able structure by using a classical 
many-body potentialE2l. This discrepancy demonstrates 
the importance of incorporating the electronic structure 
in determining cluster geometries. The lowest energy 
structure found for Ags is a compressed trigonal bipyra- 
mid, which has lower energy (Ai5=0.17 eV) than a planar 
capped rhombus. In previous studies j-tha,jdanar struc- 
ture has been obtained as ground statell3llil't3 but the en- 
ergy difference between these two isomers is rather small 
(Ai;=0.31 eV in Ref.l2 and AE= 0.003 eV in Ref.l4). It 
is noted that the experimental ESR spectra of Ags can be 
interpreted by a geometrical structure of deformed trig- 
onal bipyramid. Two isoenergitic structures, a bicapped 
tetrahedron and a pentagonal pyramid are found for Agg, 
with Aii^=0.05 eV. The bicapped tetrahedron is more sta- 
ble but this conclusion depends sensitively on the choice 
of empirical parameters. In Ref.l4, these two structures 
are also found to be very close in energy {AE — 0.06 eV) 
but the pentagonal pyramid is ground state. Accord- 
ing to the theoretical HOMO-LUMO gap of these two 
isomers, we suggest that the bicapped tetrahedron is a 
better candidate since its HOMO-LUMO gap (1.05 eV) is 
much smaller than that obtained for pentagonal pyramid 
(2.67 eV), whereas experimental gap is about 0.34 eVQ. 
The pentagonal bipyramid is obtained as lowest energy 
structure for Agy. The tricapped tetrahedron is a locally 
stable isomer with AE^ =0.48 eV, while the Ai? in Ref.l4 
for the same isomer is 0.41 eV. For silver octamer, a bi- 
capped octahedron iSrOur ground state structure, which 
is also found for CustS. This near-spherical configura- 
tion can be understood by electronic shell model. The 
closure of electronic shell at n = 8 might give rise to 
a spherical charge density distribution, which favors the 
formation of spherical atomic arrangement. The tetra- 
capped tetrahedron is predicted as metastable isomer in 
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Ref.l4 but it is unstable upon relaxation in our simula- 
tion. Another spherical-like structure, square antiprism 
(D4d) is found as a local stable isomer with Ai?=0.99 eV. 
For Agg, the ground state structure is a bicapped pen- 
tagonal bipyramid. Its energy is lower than that of the 
tricapped trigonal prism (Csv) by 0.59 eV and than that 
of capped square antiprism (C2v) by 1.01 eV. In Ref.l4, 
bicapped pentagonal bipyramid is also found as ground 
state and the energy difference IS.E for the two structural 
isomers is 0.73 eV and 0.22 eV respectively. 

Table II Comparison of structural properties (bond 
length, bond angle), average binding energies Eh/n, and 
vertical ionization potentials (IP) of smalLJVg- (n = 
3 — 9) clusters with ab initio calculationsll3~t3. The 
definition of structural parameters r, a, h, etc. for 
smaller Ag3_5 clusters is chosen according to Ref. [14]; 
the bonds for Agg-g are defined by their lengths in 
Ref. [14] in a declining sequence, a denotes our present 
tight-binding calculation. & to e are previous ab initio 
calculations based on relativistic effective core potential 
configuration (RECP): 6- modified coupled pair function 
(MCPF)E3; c-niulti*eference singles plus doubles confimi- 
ration (MRSDCI)lla; d-configuration interaction (CI)M; 
e-relativistic effective— core potential density functional 
theory (RECP-DFT)^. 



Ag3, obtuse triangle (C2v) 





r (A) 


a (deg) 


Eb/n (eV) 


IP (eV) 


a 


2.659 


66.8 


0.82 


5.65 


b 


2.709 


69.2 


0.80 


5.59 


c 


2.720 


63.7 


0.90 


5.90 


d 


2.678 


69.1 


0.86 


5.74 


e 


2.627 


70.4 


0.84 







Age, bicapped tetrahedron (C2v) 
n (A) ra (A) r-3 (A) n (A) rs (A) Eb/n{eV) 


IP(eV) 


a 
d 


2.931 2.875 2.766 2.653 2.661 1.65 
2.976 2.859 2.783 2.751 2.672 1.49 


6.71 
6.23 






Age, pentagonal pyramid (Csv) 
n (A) ra (A) Et/n{ey) 


IP (eV) 


a 
d 


2.984 2.539 1.65 
2.828 2.740 1.50 


7.92 
7.00 






Ag7, pentagonal bipyramid (Deh) 
r-i (A) ra (A) £;i,/n(eV) 


IP (eV) 


a 
d 


2.879 2.858 1.87 
2.815 2.806 1.71 


5.95 
5.91 






Ag8, bicapped octahedron (Dad) 
n (A) ra (A) rg (A) (A) Eb/n{eY) 


IP(eV) 


a 
d 


3.140 2.812 2.941 2.661 2.03 
2.973 2.812 2.804 2.735 1.80 


7.10 
6.80 




Agg, bicapped pentagonal bipyramid (C2v) 
n (A) ra (A) ra (A) (A) rg (A) 


a 
d 


2.989 2.936 2.993 2.888 
2.934 2.887 2.881 2.836 


2.856 
2.786 


a 
d 


re (A) ry (A) Eb/n{e\) IP(eV) 
2.849 2.745 2.09 5.74 
2.766 2.752 1.77 5.10 





Ag4, rhombus (Dah) 





r(A) 


a (deg) 


Eb/n (eV) 


IP (eV) 


a 


2.731 


56.6 


1.21 


6.86 


b 


2.862 


57.6 


1.11 


6.54 


c 


2.870 


55.5 


1.83 


6.40 


d 


2.800 


56.4 


1.20 


6.60 


e 


2.740 


57.2 


1.11 





Ag5, deformed trigonal bipyramid (Cav) 





r (A) 


a (deg) 


h/2 (A) 


Eb/n{eV) 


IP (eV) 


a 


2.749 


67.5 


2.34 


1.38 


5.88 


b 


2.858 


65.8 


2.39 


1.16 




d 


2.709 


67.8 


2.33 


1.28 


5.95 



Ag5, planar capped rhombus (Cav) 





ri (A) 


ra (A) 


ra (A) 


r4 (A) 


Eb/n{eV) 


IP (eV) 


a 


2.851 


2.736 


2.740 


2.668 


1.32 


6.20 


b 


2.842 


2.842 


2.842 


2.842 


1.22 


6.18 


d 


2.812 


2.801 


2.760 


2.759 


1.28 


6.20 
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From the above discussions for small silver clusters, 
we find that the agreement between TB model and ab 
initio calculations is satisfactory, particularly consider- 
ing the simplicity in the present tight-binding scheme. 
Therefore, in the next part, we shall use this model to 
study larger clusters with n > 10 for which the global 
minimization with ab initio molecular dynamics is much 
more expensive. 

B. Structures of Ag,i with 10 < n < 21 




FIG. 2. Lowest-energy structures for Ag„ {n — 10 — 21) 
clusters. See text for description of structures. 

The lowest energy structure of Ag„ (n = 10 — 21) ob- 
tained from GA-TBMD simulation is shown in Fig 2. 
The most stable structure of Agio is a deformed bicapped 
squap.antiprism (D4d), which is similar to that found for 
Cuiollj. Starting from Agn, the ground state structures 
of Ag„ clusters are based on icosahedral packing, except 
for Agi4. Many other capped polyhedral structures are 
obtained as local isomers for Ag„ with n =10-21 but it is 
not necessary to describe them herein. As shown in Fig. 2, 
the structures of Agn, Agi2 are the uncompleted icosa- 
hedron with lack of one or two atoms. An Jahn- Teller 
distorted icosahedron is then formed at Agi3. Follow- 
ing the icosahedral growth sequence, the lowest energy 
structures of Agis, Agig, Agiy is the icosahedron capped 
with 2, 3, 4 atoms respectively. The capped single icosa- 
hedron transits into an truncated double icosahedron at 
Agi8 and a complete double icosahedron at the Agig. 
Base on the double icosahedron structure, the structures 
of Ag2o and Ag2i is formed by one and two atoms capped 
on that of Agig. However, an exception occur at Agi4, 



for which we found a fcc-like packing with 4-1-4-1-4 lay- 
ered structure. This ellipsoid structure is more stable 
than a spherical capped icosahedron structure by 0.03 
eV. It is worth noted iiat 14 is a magic size predicted by 
a ellipsoid shell modeflcS. 

C. Crossover from "electronic order" to "atomic 
order" 

The concept of "particle order" (or "atomic order") 
and "wave order" (or "electronic order") in cluster 
physicg23 have been introduced in order to explain the 
magic number effect in the inert gas clusters with atomic 
shell and alkali metal clusters with electronic shell. The 
noble metal cluster is a mixture of the atomic core in- 
volving the relatively localized d electrons and the more 
delocalized s valence electrons. Therefore, it might be a 
intermediated system from these two extreme limit and 
exhibit features come from both the two orders. The 
equilibrium structures of small silver clusters (Aga— Agg) 
from our calculation are similar to those of alkali-metal 
clusters^. In contrast, the icosahedron growth sequence 
is obtained for clusters starting frnm Agn that has also 
been found for noble gas clustersH In smaller silver 
clusters, the 5s valence electrons are dominant in deter- 
mining the cluster property and the d states are signifi- 
cantly lower-lying and contribute much less to the cluster 
bonding. Therefore, these small clusters should exhibit 
certain alkali-metal-like behavior in both structural and 
electronic aspects. As the clusters size increase, the con- 
tribution of d electrons to cluster bonding become more 
and more important. The bonding energy fropi d elec- 
trons is roughly related to the d band widthcfl, which 
is approximately proportional to the square fQot of the 
average coordination number of the clustersE3. Conse- 
quently, the clusters tend to adopt the more compact 
structures such— as icosahedron, which is similar to no- 
ble gas clustersE3. The switch of structural pattern from 
alkali-metal like to noble gas like at around n=10 can be 
identified as a crossover from "electronic order" towards 
"atomic order" in the silver clusters. This alternation is 
related to the overlap of 4d and 5s electronic states which 
we will discussed latter. However, our further study show 
that the shell structure of s electrons still dominates the 
electronic properties such as IPs, HOMO-LUMO gaps 
of the silver clusters although the geometrical structures 
has taken "atomic order" . We argue that the "electronic 
order" and "atomic order" can coexist in coinage metal 
clusters. 



D. Size dependence of relative stabilities 

The second differences of cluster binding energies de- 
fined by A2E{n) = Eb{n + 1) + Eb{n - 1) - 2Eb{n) is 
calculated and plotted in Fig. 3. In cluster physics, it 
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is well known that the A2E{n) is a sensitive quantity 
which reflects the stability of clusters and can be directly 
compared to the experimental relative abundance. Three 
major characteristics can be found in the Fig. 3: (i) even- 
odd alternation of A2-E(n) with n = 2 — 6, 15 — 21; (ii) 
particular high peak at Agg, Agig; (iii) other maxima 
at odd size like Agia and Agn. The first effect can be 
related to the even-odd oscillation of HOMO energy and 
HOMO-LUMO gap in silver clusters, which is due to elec- 
tron pairing effect. The articular stable clusters such as 
Agg, Agig corresponds to the magic number in electronic 
shell model. However, the even-odd oscillation in A2E{n) 
from Agio to Agi4 and the maximum at magic size Ag2o 
have not been observed in our calculation. In stead, some 
odd-sized cluster as Agn, Agia become maxima in Fig. 3. 
These phenomena can be attributed to the geometrical 
effect. The closing of geometrical shell of icosahedron 
at Agi3 will enhance the stability of such clusters and 
reduce the relative stability of their neighboring clusters. 



experimental probe of jgeometrical and electronic struc- 
ture of copper clusterstj. They found both jellium-like 
electronic behavior and icosahedral geometrical structure 
in copper clusters. In a expecipiental studies of mass 
spectra of ionized silver clustergj, dramatic even-odd os- 
cillation as well as substantial discontinuities at electronic 
magic number 8, 20 (n = 9, 21 for cationic clusters) are 
found. The discrepancy between present theoretical re- 
sult and experiment may be partially attributed to the 
effect of ionization on the cluster stability. Since the ex- 
perimental mass spectra distribution is recorded for ion- 
ized clusters Ag+ , it is possible that the charging on the 
cluster can significantly alter|-|the geometrical and elec- 
tronic structure of the clusterElllJ. 



IV. ELECTRONIC PROPERTIES VS CLUSTER 
SIZE 

A. Size evolution of electronic band 
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FIG. 3. Second differences of cluster binding energies 
AE{n) = [i5i,(Ag„_i) + £5(Ag„+i)] - 2Et{Ag„) as a func- 
tion of cluster size n for n = 2 — 21. Both electronic shell 
effect at n = 2, 8, 18 and geometrical shell effect at Agis can 
be identified. See text for details. 

The simultaneous appearance of those three features 
in the A2E{n) demonstrates that the structure and sta- 
bility of a silver cluster is determined by both electronic 
structure and atomic configuration. Either electronic or 
geometrical effect is enhanced if the corresponding shell 
structure is completed. This argument is supported by a 



We investigated the cluster electronic properties via 
the electronic density of states (DOS). In Fig. 4, we 
present the total spd electronic DOS for Ag2, Agg, Agi3 
along with bulk DOS of fee crystal silver from TB cal- 
culation in reciprocal space. Generally, the total DOS 
is composed by the relatively compact d states and the 
more expanded sp states. In smallest clusters such as 
Ag2, the d and sp bands are clearly separated. The sp 
states shows discrete peaks originated from symmetrical 
splitting of atomic orbital levels, while the d band is low- 
lying and considerably narrower than the bulk d band. In 
contrast to even-odd behavior and shell occupation of s 
electrons, the evolution of d states from smallest clusters 
towards bulk solid is a monotonic broaden of band width. 
As the cluster size increases, both d and sp levels grad- 
ually broaden, shift, overlap with each other, and finally 
come into being bulk electronic band. The DOS of Agg 
still has molecular-like some discrete sp peaks but these 
electronic spectra peaks tend to overlap and form con- 
tinuous band. In Agia, the sp states have developed into 
several subbands and the d band has overlapped with sp 
states. Although the DOS of Agi3 is roughly similar to 
the bulk limit, the width of d band is still considerably 
narrower than the bulk d band width and the fine struc- 
ture of sp electronic spectra is somewhat different from 
bulk sp band. This fact suggests that the bulk-like be- 
havior emerge at around Agia. We have also studied the 
electronic states of Agss with icosahedral and cubocta- 
hedral structures by using present tight-binding scheme 
with local minimization. The DOS for both of them are 
much closer to bulk band. In a previous experimental 
study of jphotoelectron spectra of silver clusters up to 
60 atomsB, it is found that the ultraviolet photoelectron 
spectroscopy (UPS) of smallest Ag„, i.e., 2 < n < 10 is 
different from bulk UPS and changes sensitively on clus- 
ter size. The size dependent variation of UPS for Ag„ 
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with n < 10 becomes more gradual and the UPS of Aggo 
is aheady very close to that of solid silver. 



B. HOMO-LUMO gap 
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FIG. 4. Density of states (DOS) of Ag„ (n = 2, 8, 13) 
clusters vs cluster as well as bulk DOS in fee crystalline (with 
Gaussian broaden of 0.02 eV). 

To further clarify the size evolution of the overlap of 
d and sp bands in the small silver clusters, we have ex- 
amined the energy separation Asd between the highest 
molecular orbitals belong to d states and lowest molecu- 
lar orbitals from s states. The calculated Agd decrease 
rapidly from 5.20 eV for Ag2, to 1.69 eV for Ags and 
then to 0.10 eV for Agg. Finally, the d and sp band 
merge in Ag„ clusters with n > 9. The overlap between s 
and d band is believed to induce the icosahedral growth 
sequence starting from Agn and weaken the even-odd 
oscillation in HOMO-LUMO gaps and IPs of Ag„ clus- 
ters with ri > 10. However, the overlap between sp and 
d states is small and the cluster HOMO is still located 
in the s-like electronic states. Consequently, the contri- 
bution from d electrons in Ag„ with n > 9 shows more 
importance to cluster bonding although the electronic be- 
havior close to Fermi energy such as HOMO-LUMO gap, 
ionization potentials is still dominated by s orbitals. 



FIG. 5. HOMO-LUMO gap A (eV) vs cluster size n. The 
theoretical values are labeled by solid square connected with 
solid line and the experimental values in Ref.[7] are labeled 
by open circles. Electronic shell effect for n — 2, 8, 18 can be 
clearly identified. 

An important electronic property of a cluster is the 
gap between highest occupied molecular orbital (HOMO) 
and lowest unoccupied molecular orbital (LUMO). In the 
case of magic cluster, the closure of electronic shell shall 
manifest itself in particularly large HOMO-LUMO gap. 
This effect was demonstrated experimentally for small 
even-sized silver and, cppper clustersQ and theoretically 
for copper clusterMN. The theoretical HOMO-LUMO 
gap of Ag„ (n = 2 — 21) along with experimental gap 
of small clusters Ag„ (n = 2, 4, 6, 8, 10)0 are shown 
in Fig. 5. Even-odd oscillation up to Agig as well as 
the particularly large HOMO-LUMO gap at Ag2, Ags, 
and Agi8 are obtained. As compared to the experimen- 
tal results for small Ag„ with even size, the present TB 
calculation has systematically overestimated the HOMO- 
LUMO electronic gap by about 0.5 eV. But the size de- 
pendent variation of experimental gaps and magic effect 
in HOMO-LUMO gaps at n =2, 8 are well reproduced. 
The even-odd alternation for n > 16 and magic effect 
of Ag2o have not been obtained in our calculation. We 
suggest these are probably due to the geometrical effect, 
since the HOMO-LUMO gap of cluster depends sensi- 
tively on cluster structureEH. In a previous study of 
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HOMO-LUMO gaps of copper cluster£3, the maxima 
gap at Ag8 and Ag2o is found but the even-odd alterna- 
tion of electronic gap and magic effect for Agig have not 
been obtained. 



C. Ionization potential 

The vertical ionization potentials (IPs) of clusters are 
evaluated from the highest occupied molecular orbital 
(HOMO) energy of neutral clusters according to Koop- 
man's theorem. In Fig. 6, the calculated IPs of Ag„ 
up to n = 21 is compared with the IP values mea- 
sured by C.JackschathB, the prediction by metallic spher- 
ical droplet models, and the size dependent HOMO 
level (in arbitrary units) of alkali-like W^tal clusters by 
Clemenger-Nilsson ellipsoid shell modelElO. In compari- 
son with experimental values in Ref.[5], the present TB 
calculation has almost reproduced the size dependence 
of IPs for silver clusters up to 21 atoms except that the- 
oretical calculation has overestimated the IP values of 
some magic clusters such as Ag2, Agg, and Agig. Two 
important size dependent behaviors are found in Fig. 6; 
(i) dramatic even-odd alternations where clusters with 
even number of s valence electrons have higher IPs than 
their immediate neighbors; (ii) particular higher IP val- 
ues at the magic clusters such as Ag2, Agg, Agig, Ag2o 
and probably Agi4. The even-odd variations can be at- 
tributed to electron pairing effect. Odd(even)-sized clus- 
ters have an odd (even) total number of s valence elec- 
trons and the HOMO is singly(doubly) occupied. The 
electron in a doubly occupied HOMO feels a stronger 
effective core potential since the electron screening is 
weaker for the electrons in the same orbital than for in- 
ner shell electrons. Therefore, the binding energy of a 
valence electron in a cluster of even size cluster is larger 
than that of odd one. It is also interesting to note that 
the size dependence of IPs from TB model is almost in 
full accordance to Clemenger-Nilsson shell model. The 
magic effect at Ag2, Agg, AgM, Agig, and Ag2o pre- 
dicted by electronic shell model is well reproduced even 
though the cluster geometries and s-d interaction has 
been considered. On the other hand, one can found that 
both the theoretical and experimental IP values for silver 
clusters can be roughly described by a classical electro- 
static model which take the cluster as a metallic spherical 
droplet. This fact further suggests that the electronic be- 
havior of silver clusters are predominantly s-like and the 
cluster can be approximated to be a structureless jellium 
sphere with shell-like electronic levels. 
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FIG. 6. Vertical ionization potentials (IPs) vs cluster size 
n. The measured IP values reported in Ref. [5] are labeled by 
open circles connected with dashed line, the theoretical values 
from TB calculation are labeled by open squares connected 
with solid line. Electronic shell closing exhibit as pronounced 
drops of IP at n = 2, 8, 14, 18, 20. See text for details. 

V. CONCLUSION 

In conclusion, we have shown that the present TB 
model is able to describe the geometrical and electronic 
structures of silver clusters. By using a combined GA- 
TBMD technique, the lowest energy structures, binding 
energies, electronic states, s-d separation, HOMO-LUMO 
gap, vertical ionization potentials are obtained and com- 
pared with experiments. The main theoretical results can 
be summarized in the following points: 

(1) The structures of small silver clusters is determined 
by s electrons and similar to those of alkali-metal clusters. 
The contribution of d electrons to cluster bonding become 
more important as cluster size exceed 10 atoms. The 
icosahedral growth pattern starts in the Ag„ with n > 11, 
which can be identified as "atomic order" clusters like 
noble gas clusters. 

(2) The electronic and geometrical shell structure coex- 
ist in silver clusters and take effect in the clusters simul- 
taneously. The electronic shell effect on cluster electronic 
properties has been found by present TB calculation, in 
which the effect of geometrical structures and d electrons 
are directly incorporated. The silver clusters with closed 
electronic shell {n —2, 8, 14, 18, 20) show more pro- 
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nounced electronic characteristics while the geometrical 
effect is enhanced as the icosahedral shell completes at 

Agl3. 

(3) Due to the pair occupation of s valence electrons 
on molecular orbitals, silver clusters show even-odd al- 
ternation in their relative stability, HOMO-LUMO gap, 
ionization potential. However, the even-odd effects can 
be disturbed by the sd overlap or the geometrical struc- 
tures of clusters. 

(4) The density of electronic states of smaller silver 
cluster, e.g., n < 10, is composed by discrete sp expanded 
band and a narrow d band. The bulk-like feature in DOS 
start at around Agi3 and the bulk limit can be roughly 
reached by n = 55. 

The present study shows that both the geometrical and 
electronic effects should be considered in order to achieve 
a complete description of coinage clusters. Therefore, ab 
initio molecular dynamics or TBMD are essential to elu- 
cidate the interplay between geometrical and electronic 
structures of these clusters. The further works should 
include the larger clusters and extend the TB model to 
other transition metal elements. 
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